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RESUMEN 

Presentamos modelos de evolution quimica del disco Galactico con diferentes 
rcndimicntos dcpcndicntes de Z. Encontramos que una tasa moderada de 
perdida de masa en estrellas masivas de metalicidad solar produce un excelente 
ajuste con los gradientes de C/H y C/O del disco de la Galaxia. El mcjor 
modelo reproduce: las abundancias de H, He, C, y O derivadas de hneas de 
recombination en Ml 7, las abundancias protosolares y las relaciones C/O- 
O/H, C/Fe-Fe/H y O/Fe-Fe/H derivadas de estrellas de la vecindad solar. El 
acuerdo del modelo con las abundancias protosolares implica que el Sol se 
origino a una distancia galactocentrica similar a la actual. El modelo para 
r = 3 kpc implica que una fraction de las estrellas en la direction del bulbo 
se formo en el disco interno. Nuestro modelo reproduce la relation C/O-O/H 
dcrivada de regiones H n extragalacticas en galaxias espirales. 

ABSTRACT 

We present chemical evolution models of the Galactic disk with different .In- 
dependent yields. We find that a moderate mass loss rate for massive stars 
of solar metallicity produces an excellent fit to the observed C/H and C/O 
gradients of the Galactic disk. The best model also fits: the H, He, C, and 
O abundances derived from recombination lines of Ml 7, the protosolar abun- 
dances, and the C/O-O/H, C/Fe-Fe/H, and O/Fe-Fe/H relations derived from 
solar vicinity stars. The agreement of the model with the protosolar abun- 
dances implies that the Sun originated at a galactocentric distance similar 
to the one it has. Our model for r = 3 kpc implies that a fraction of the 
stars in the direction of the bulge formed in the inner disc. We obtain a good 
agreement between our model and the C/O versus O/H relationship derived 
from extragalactic H n regions in spiral galaxies. 

Key Words: galaxies: abundances — galaxy: bulge — galaxies: evolution — 
H II regions: M17, Orion nebula — ISM: abundances — Sun: abundances 



1. INTRODUCTION 

The comparison of detailed Galactic chemical evolution models, GCE mod- 
els, with accurate abundance determinations of stars and gaseous nebulae pro- 
vides a powerful tool to test the chemical evolution models and the accuracy 
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of observational abundance determinations of stars of different ages and of 
H II regions located at different galactoccntric distances. 

In this paper we will compare our models with stellar and H n regions 
abundances to test if the H n region abundances derived from recombina- 
tion lines agree with the stellar abundances, in particular with the protosolar 
abundances that correspond to those present in the interstellar medium 4.5 
Gyr ago. Also our GCE models can be used to constrain the C yields for 
massive stars, the C yield is not well known and we will vary it to obtain the 
best fit between our GCE models and the observational data. 

Carigi and Peimbert (2008, hereinafter Paper I) presented chemical evo- 
lution models of the Galactic disk for two sets of stellar yields that provided 
good fits to: a) the O/H and C/H gradients (slope and absolute value) de- 
rived from H II regions based on re combination lines (E steban et al. 2005) 
and including the dust contribution (lEsteban et al 1 ll998l) . and b) the AY/AZ 
value derived from the Galactic H II region M17, and the primordial helium 
abundance, Y p obtained from metal poor extragalactic H n regions (Peimbert 
et al. 2007). In Paper I, based on our GCE models and combined with the 
constraints available, we were not able to discriminate between the stellar evo- 
lution models assuming high wind yields for massive stars, HWY, and those 
assuming low wind yields for massive stars, LWY. 

Previous works have focused on the test of stellar yields using GCE models 
constrained by chemical gradients obtained by different methods, gradients 
that in general show similar slopes but a c o nsiderable sp r ead in the abso 



lute O/H ratios, e.g. iPrantzos et al.l (|1994f ); ICarigil (119961 ); IChiappini et al 



(|2003aj ): iRomano et al.1 (|2010l ). To test the stellar yields it is necessary to 
have good absolute abundance values of stars and H n regions. To deter- 
mine the O/H abundances of H n regions in irregular and spiral galaxies 
many methods have been used in the literature. Most of them have been 
based on fitting photoionizcd models to observations or by determining the 
electron temperature, T, from the 4363/5007 [O III] ratio directly from obser- 
va tions. A comparison of m any of the different methods used has been made 
by Kewlev fc Ellison! (|2008f ). They find that the O/H differences derived by 
different methods between two given H II regions amount to 0.10 — 0.15 dex. 
Alternatively for all the methods the absolute difference for a given H II re- 
gion is considerably larger reaching values of 0.7 dex for extreme cases (see 
Figure 2 in their paper and the associated discussion) . Most of the differences 
among the various calibrations are due to the temperature distribution inside 
the nebulae. In this paper we will use only abundances of H II regions based 
on recombination lines of H, He, C, and O, these lines depend weakly on the 
electron temperature, they are roughly proportional to 1/T, therefore the rel- 
ative abundances among these four elements are practically independent of 
the electron temperature. 

There are two frequently used methods to derive C and O gaseous abun- 
dances from H II regions: a) the most popular one based on collisionally 
excited lines (or forbidden lines) and the T(4363/5007) temperatures, the FL 
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method, and b) the one based on C and O recombination lines, the RL method. 
The RL method produces gaseous O and C abundances higher by about 0.15 
to 0.35 dex than the FL method. The RL method is almost independent of 
the electron temperature, while the FL method is strongly dependent on the 
electron temperature. It is possible to increase the FL abundances under the 
assumption of temperature inhomogeneities to reach agreement with the RL 
values. The temperature distribution can be characterized by the av erage tem- 
perature, Th, and the mean square temperature variation, t 2 , (e.g. iPeimbert 



19671 ; IPeimbert et al.l |2002| ) . The t 2 values needed to reach agreement be- 



tween the RL and the FL abundances are in the 0.02 to 0.05 range, while the 
photoionization models predict typically t 2 values in the 0.003 to 0.01 ra nge, 
this discrepancy needs to be sorted out (e.g. IPeimbert fc Peimbert II 2 1 ll . and 
references therein) . 

Paper I is controversial because the C /H and O /H gaseous abundances of 
the H II regions have been derived from recombination lines (that is equivalent 
to the use of t 2 ^ 0.000 and forbidden C and O lines) and the assumption 
that 20% of the O atoms and 25 % of the C atoms are trapped in dust grains 
(0.08 dex and 0.10 dex respectively). These assumptions increase the O/H 
ratio by about 0.25 to 0.45 dex relative to the gaseous abundances derived 
from T(4363/5007), the forbidden O and C lines, and the assumption that 
t 2 = 0.00. 

Due to the controversial nature of the H II region abundances used in 
Paper I and that we were not able to discriminate between the two sets of 
stellar yields adopted we decided to test our GCE models further by includ- 
ing additional observational constrains: a) the Asplund et al. (2009) proto- 
solar abundances that provide us with the O/H, C/H, Fe/H, and AY/AO 
in the interstellar mediu m 4.5 Gyr ago, b) the C/H, O/H, and Fe/H by 



Bensbv &: Feltzina (|2006l ) for young F and G stars of the solar vicinity, c ) 



the O/H, C/H, and He/H derived from B stars by iPrzvbilla et al] (|2008f ). 
and d) throughout this paper for all the H n regions we will use abundances 
derived from recombination lines and to obtain the total abundances we will 
increase the gaseous abundances by 0.10 dex in C and 0.12 dex in O to take 
into account the fraction of atoms trapped in dust grains, with the exception of 
the me tal poor irregular ga l axies for which we will use an 0.10 dex depletion 



for O (Esteban et al 



1998; 



Mesa-Delgado et al 



2009; 



Peimbert fc Peimbert 



2010). 



We also decide d to compare our bes t mod els with the C/O versus O/H 
results derived by Esteban et al. ( 2002 . 2009h from bright H II regions in 
nearby spiral galaxies based on recombination lines and to make a preliminary 
discussion of a comparison between our models f or r = 3kpc and the stars 
in the direction of th e galactic bulge obtained by iBensbv et al.l (|2010al ) and 



Zoccali et al. 



(120081) . 

The symbols C, O, X, Y, and Z represent carbon, oxygen, hydrogen, 
helium, and heavy element abundances by unit mass respectively; while C/H, 
O/H, Fe/H, C/O, C/Fe, and O/Fe represent the abundance ratios by number. 
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In Section [2] we discuss the general properties of the chemical evolution 
models, we discuss infall models for the Galaxy with two sets of stellar yields, 
the HWY and the LWY. In Section [3] we show the prediction of the current 
abundance gradients for the interstellar medium (ISM) of the Galactic disk 
and compare them with Galactic H n region abundances derived from re- 
combination lines that include the dust correction, in addition for the solar 
vicinity we present the chemical history of the ISM and compare it with the 
chemical abundances of stars of different ages; we define the solar vicinity as 
a cylinder perpendicular to the galactic plane, centered in the Sun, with a ra- 
dius of 0.5 kpc, that extends into the halo to include the stars in the cylinder. 
In Section [4] we compare our chemical evolution models with the protosolar 
chemical abundances and with those of the H n region M17. Based on the 
comparison between the observations and the models in Section [5] we present 
a new Galactic chemical evolution model, with intermediate mass loss due to 
interstellar winds, IWY, that produces considerably better adjustments with 
the observations. In Section [6] we compare the IWY model with additional 
observations, those provided by extragalactic H n regions in spiral galaxies, 
and those provided by stars in the direction of the Galactic bulge. The con- 
clusions are presented in Section [7] A preliminary account of some of th e 



results included in this paper was presented elsewhere iPeimbert et al.l (|2010D 



2. CHEMICAL EVOLUTION MODELS WITH HIGH WIND YIELDS 
AND LOW WIND YIELDS 

We present chemical evolution models for the Galactic disk. The models 
have been built to reproduce the present gas mass distribution, (see Fig. 1 , left 
panel) and the present-day O/H values for H II regi ons in the Galaxy for 6 < 
r(kpc )< 11, (see Fig. 1, upper right panel) listed by Garria-Roias k, Estebanl 



(|2007h . In what follows we describe in detail the characteristics of the models. 

i) In the models the halo and the disk are projected onto a single disk 
component of negligible width and with azimuthal symmetry, therefore all 
functions depend only on the galactocentric distance r and time t. 

ii) The models focus on r > 4 kpc, because the physical processes associ- 
ated with the Galactic bar are not considered. 

iii) The models a re based on th e standard chemical evolution equ ations 



origin a lly written by iTinslevi ( 19801) and wide ly used to date. See e.g. iPagel 
(j2009h . lMatteuccil |2000h . and IPrantzol |2008h . 



iv) The age of the models is 13 Gyr, the time elapsed since the beginning 
of the formation of the Galaxy. 

v) The models are built based on an inside-out scenario with infalls that 
assume primordial abundances (Y p = 0.2477, Z = 0.00, Peim bert et al. 2007) 



The a dopted double infall rate is similar to that presented bv lChiappini et al 



( 19971 ) . as a function of r and t, and is given by IR(r,t) = j4(r)e - '/ Thal ° + 
B(r)e~( t ~ 1Gyr ^ Tdi!3k , where the halo formed in the first Gyr with a timescale, 
Thaio == 0.5 Gyr, and the disk started forming immediately after with longer 
timescales that depend on r, rjisk = (8 x r/r — 2) Gyr. We adopt 8 kpc 
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for the galactocentric distance of the solar vicinity, r . The variables A(r) = 

lOMopc" 2 X e -(r(kpc)-r e )/3.5kpc and _ ^MqPC" 2 X e -(r(kpc)-r Q )/3.5kpc 

are chosen to match the present-day mass density of the halo and disk com- 
ponents in the solar vicinity that amount to 0.5 and 49.5 TsAqpc~ 2 , respec- 
tively, where the mass lost by the stars and the halo gas have been incor- 
porated into the disk. Moreover with the A(r) and B(r) variables the mod- 
els reproduce the radial profile of the total mass in the Galaxy, M tot (?') = 



5OM pc" 



x e 



— (r(kpc)- 



re)/3-5k P c (jFenner fc Gibson! I2003T ). 



vi) The models assume the Initial Mass Function (IMF) proposed by 
Kroupa et al.l (119931 ). in the mass interval given by 0.08 — 80.0 M Q mass range 
for Z > 10 -5 , and in the 9.0-80.0 M Q mass range for Z < 10~ 5 . We consider 
that Population III do not include objects with less t han 9.0 and that 



the change from Pop III to Pop II. 5 occurs at 10 Zq (jAkerman et al.ll2004 ; 



Bromm fc Lars on 2004). 



vii) The models include a star formation rate that depends on time and 



gala ctocentric distanc e , SFR (r, t) 



,t) (M nn . s +M stars ) 0A (r, t), taken 



from lMatteucci et al.l (I1989T ) and iMatteucci fc Chiappinil (119991 ) . where v is a 
constant in time and space. This SFR formula considers the feedback between 
gas and stars. In our models v is chosen in order to reproduce the present-day 
radial distribution of gas surface mass density. We adopted v values of 0.019 
and 0.013 for the HWY and the LWY models respectively. We assumed a v 
value 5 times higher during the halo formation than that adopted for the disk. 
These v values combined with the infall rate adopted reproduce the chemical 
abundances shown by halo and disk stars (see Fig. 2). 

viii) The only difference between the HWY and the LWY sets is the as- 
sumed mas s-loss rate due to st e llar w inds by massive stars with Z = 0.02, see 
Figure 1 of ICarigi fc Peimbertl (|2008l ). All stellar yields are metal dependent. 
We interpolate and extrapolate linearly the yields by mass and metallicity. 

The HWY set includes: 

A) For ma ssive stars (MS), those with 8 < m/M Q < 80, the yields by: a) 
Hirschil (|2007l ) for Z = 10~ 8 (wit h rotation velocity betwee n 500 and 800 km/s, 
depending on stellar mass); b) iMevnet fc Maederl (120021) for Z = 10~ 5 and 
Z = 0.004 (with rotation velocity = 300 km/s): cl lMaederl (|l992h for Z = 0.02 
(high mass- loss rate with no rotation velocity); d) Since Fe is a chemical 
element used as observational constraint in most of the chemical evolution 
model s and Fe yield is no t compute d bv iHirschil (120071) . Mevnet fc Maeder 

1995J) only for the 



Weaver 



(|2002l ). and iMaederl (|1992l ). we adopt IWooslev 
Fe yields (Models B, for 12 to 30 M Q ; Models C, for 35 to 40 M Q ; while for 
m > 4OM , we extrapolated the m — 40M Q Fe yields, following Carigi & 
He rnandez 2008 ). Our mo dels, that include the Fe yields and the O yields 
bv lMevnet fc Maederl (|2002l ). can reproduce the O/Fe-Fe/H trend of the halo 
stars (see Fig 2). 

B) For low and intermediat e mass stars (LIMS) , those with 0.8 < to/M w < 



8, we have used the yields by iMarigo et al.l (119961 Il998l ) and iPortinari et al 
|l998h for the Z = 0.004 to Z = 0.02 range. 
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Thielemann et al 



(11993) 

19831 ). A fraction, Abin, 



C) For Type la SNe we hav e used the yields by 
in the the SNIa formulation of iGreggio fc Renzini 
of the binary stars of LIMS with a total mass between 3 and 16 M Q are 
progenitors of SNIa. We assumed Abin = 0.066 for HWY model, while Abin = 
0.094 for LWY models. These fractions are needed to reproduce the age- 
[Fe/H] relation of the disk stars, the O/Fe-Fe/H relation of the disk stars of 
solar vicinity, and the protosolar Fe/H value for each model (see Fig. 2). 

In the LWY set we have up dated the yie l ds of massive stars only for 
Z = 0.02 assuming the yields bv iHirschi et al.1 (|2005l ) with rotation velocity 
= 300 km/s. The rest of the stellar yields are those included in the high wind 
set. 



ix) For MS the models inclu de the stellar l i fetime s given bv lHirschil ([20071 ) . 



Mevnet fc Maederl ([20021 ). and IHirschi et al.1 (120051). w hile for LIMS include 



the main sequence lifetimes from ISchaller et al.l (|1992[) . 

x) The models do not consider any type of outflows from the Galaxy due 
to its deep potential well. In addition, we discard outflows of C and O rich 
material from SNe , due t o the small C/O ratios observed in H II regions. See 



Carigi et al.l (|l995l Il999h . 



xi) We do not include radial flows of gas or stars. 

Since the solar vicinity and the Galactic disk contain stars and H II regions 
of a broad range of metallicities, our galaxy is a proper laboratory to study the 
AY/ AO and AC/ AO behavior at high Z values and to test obscrvationally 
the predictions of the HWY and the LWY models. Furthermore the main 
elements affected by the winds of massive stars are He, C, and O, therefore 
a careful comparison of the abundances predicted by the chemical evolution 
models with the observed values will permit to constrain the mass loss rate 
and consequently the yields of massive stars. 



3. THE GALACTIC H II REGIONS GRADIENTS AND THE CHEMICAL 
HISTORY OF THE SOLAR VICINITY 



In Figure 1 we show the present-day radial distribution of the gaseous 
mass and the O/H, C/H, and C/O gradients in the Galactic disk predicted 
by the HWY and LWY models. 

For Figure 1 we have chosen observational constraints that represent the 
current gas mass and the abundances of the Galactic disk. We have taken 
the average s urface density as a funct ion of the galactic radius derived from 
H I shown by Kalberla fc Kerp! ( 2009 ) in their Figure 5. We have assumed 
the abundance ratios of the gas component determined f r om H n regions 
based on recombination lines by iGarcia-Roias fc Estebanl ([20071 ) corrected 
by dust depletion. It s h ould be noted that the the H II regions used by 



Garcia-Roias fc Estebanl (|2007l ) are high density young objects that do not 
contain WR stars, consequently they have not been polluted by the evolution 
of their ionizing stars and their O/H and C/H values are representative of the 
present value of the galactic interstellar medium. 
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Fig. 1. Present-day radial distribution of: gas surface mass density (left panel), 
and ISM abundance ratios (right panels). Predictions of our chemical evolution 
models for the Galactic disk at the present time: HWY (continuous lines), LWY 
(dashed lines). Observation al data: area enc l osed b y dotted lines: average gas sur- 
face density distribution b y Kalberla fc Kerp ( 20091 ) . filled circles: H II regions, gas 
(jGarcia-Roias fc Estebanll2007l ) plus dust values (see text). 



Based on Figure 1, it can be noted that both models successfully reproduce 
the current radial distribution of gas surface mass density and the C /H and 
O/H gradients at the one a error but neither of them reproduces the C/O 
gradient for all Galactocentric distances. The HWY model adjusts the C/O 
values of H n regions for r < 7.3 kpc, while the LWY model does so for 
7.5 < r(kpc)< 9. 

In the literature there are chemical abundance determinations from recom- 
bination lines of Galactic H II regions only for 6 < r(kpc)<ll. At present, 
the evidence for the flat tening of the Galactic gradients is not conclusive: 



Vflchez fc Estebanl (|1996l ) based on H II regions found a flat O /H gradient for 



r > 14 kp c, while more re cent work based on Cephcids did not confirm this 
result, see Mattsson ( 201fj|) . Additional observations of high quality obtained 
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with the same method are needed to established the behavior of the gradients 
for r > 11 kpc. With an inside-out scenario without dynamical e ffects it is 
difficu lt to reproduce a possi ble flattening of the O/H gradient, see 



Mattsson 



(]2010l ). ICescutti et al.l (|2007f ) in an inside-out scenario reproduce the flatten- 
ing of radial gradients shown by Cepheids for r > 7 kpc assuming a constant 
surface density for the halo. In our models it is equivalent to adopt a constant 
A(r), which means that the volumetric density of the halo increases towards 
the outer parts of the Galaxy, a densi ty di stribution that we consider un - 
likely. According to iRoskar et al.l (|2008l ) and ISanchez-Blazquez et al.l (|2009| ) , 
the flattening is a consequence of stellar migration and of a break in the star 
formation rate at large radii. A third possibility is that the formation of the 
Galaxy had a modest stellar formation previous to the inside-out scenario, 
support for this idea comes from the increase of th e average age of the stars 
at large galactocentric distances (jVlaiic et al]|201ll and references therein). 

In Figure 2 we show the evolution of C/O-O/H, C/Fc-Fc/H, O/Fe-Fe/H, 
and time-Fe/H, relations, predicted by the HWY and LWY models for r = 8 
kpc. The time vs Fe/H plot is an equivalent representation to the known 
age-Z relation of the solar vicinity. 

For Figure 2 we have chosen observational constraints that represent the 
chemical history of the solar vic i nity. We have taken dwarf stars of the 
Galactic halo by lAkerman et al.l (]2004f ). as representative of the first Gyr 
of the evolution, a nd dw arf stars of the Galactic thick and thin disk by 



Bensbv fe Feltzind (|2006l ). as representative of the last 12 Gyrs of the evo- 



lution. Also in Figure 2 we present the mean a ges and Fe/H values of d isk 
stars with better determined ages presented by iNordstrom et al.l (|2004j ) in 
the lower part of their Figure 28. M oreover, we have assumed the protosolar 

abundances bv lAsplund et al.l (|2009l ). as representative of t = 8.5 Gyr. 

W e converted the chemical abundances determined by iBensbv fe Feltzing 



(|2006 ) to abundance ratios by number assuming their own solar ab undances, 



see th eir paper for references. Also, we normalized the stellar ages bv lNordstrom et al 
|2004h to the age of the model (13.0 Gyr). 

Based on Figure 2, it can be noted that both models produce a reasonable 
fit to the C/O-O/H, C/Fe-Fe/H, and O/Fe-Fe/H trends in the solar vicinity. 
From the C/O-O/H and C/Fe-Fe/H relations it can be seen that the HWY 
model predicts more C than observed in metal rich disk stars while the LWY 
model predicts less than the HWY. Alternatively both models predict a C /Fe 
plateau for Fe/H higher than solar, while metal rich stars of the thin disk show 
a C/Fe decrease. This observed trend could be explained if massive stars of 
Z > Zq are less efficient producing C than massive stars of solar metallicity. 
Stellar yields of MS and LIMS for Z > Z Pl ar e needed to have a complete 
pic ture of the evo l ution at hi gh Z (Carigill2008 ). 

Cescutti et all |2009h and iRomano et al.l |2010h focused on the C/O-O/H 



relation shown by thin disk stars of the solar vicinity and found that the C/O 
rise at high O/H values can b e explained partially with metallicity- depend 
stellar winds in massive stars ( Maeder 1992; iMevnet fe Maederl 120021) . This 
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Fig. 2. Chemical evolution models for the solar vicinity (r = 8 kpc): HWY 
(continuous lines), LWY (dashed lines). The left panel shows the C/O evolution 
with O/H. The right panels show the evolution of C/Fe and O/Fe with Fe/H, and 



the F e/H-time relation. Filled red squares: halo dwarf stars from lAkerman et al 
( 2004|V_ Filled blue and e mpty green triangles: thick-disk and thin-d isk dwarf-stars 
fromlBensbv fc Feltzin3 (|2006h . Dotted lines: protosolar values from lAsplund et al 
(2009). Open ci rcles: mean ages and Fe /H values of disk stars with better deter- 



mined ages from iNordstrom et al.l (|2004r ). Horizontal bars: internal dispersions in 
Fe/H shown by the sample. Vertical bars: age average errors. 



result is in agreem ent with our previous conclusions ( Carigi 1994 L 



1996 



good 



Carigi et al.ll2005l and Paper I), and with those bv lPrantzos et al.l (|1994[ ). 



The fit of our models to the time- Fe/H relation shown by disk stars of 
the solar vicinity is reasonable for 4 < i(Gyr)< 13, but our models cannot 
reproduce the mean behavior of older stars of the Galactic disk. Probably 
the reason is that some of these stars originated closer to the center of the 
Galaxy and migrated outwards, or belonged to satellite galaxies with different 
chemical histories that were captured by our galaxy. Most of the stars of the 
sample have ages lower than 9 Gyr, corresponding to t > 4 Gyr, and the 
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internal dispersion in Fe/H shown by old disk stars is higher than that by 
young stars. Both models in the t < 1 Gyr range adjust the Fc/H values of 
the halo stars with ages between 12 and 13 Gyr. 

The main difference between the HWY and the LWY sets is due to the 
stellar yields assumed for massive stars at Z = 0.02. The HWY assu me a rel- 



ative ly high mass-loss rate for massive stars with Z = 0.02 (yields by iMaeder 



1992J) , while the LWY as sume a relatively lo w mass- loss rate for massive stars 



with Z = 0.02 (yields bv lHirschi et al.ll2005l ). Since the mass loss rate is pro 



portional to the stellar metallicity, the efficiency of this rate i ncreases with 



metallicity and becomes important at Z ~ Zq . According to iHirschi et al 



(|2005l ) mass loss rates are a key ingredient for the yields of mass i ve sta rs and 



the rates assumed by them arc 2-3 smaller than those by Maeder ( 1992h . This 



difference between a high and a low mass-loss rate produces opposite differ- 
ences in the O and O yields, these differences occur in the pre-SN and in the 
SN phases. In the pre-SN phase massive stars are able to process He into 
C, and their stellar winds take away a lot of new C. Since the O production 
happens deeper than the C production, the wind contribution to the O yield 
is much smaller than to the C yield; consequently C constitutes the largest 
fraction of heavy elements ejected during the wind phases. While, in the SN 
phase the contribution to the C and O total yields depends mainly on the mass 
of the CO core, a small fraction of C in the CO core remains unmodified and 
is ejected during the explosion. Alternatively the O production by the SNe is 
proportional to an important fraction of the mass of this core; consequently O 
constitutes the largest fraction of heavy elements ejected during the SN stage. 
Therefore, when the initial stellar metallicity is higher, the mass loss rate is 
higher, the C yield is higher, the CO core is less massive, and consequently 
the O yield is smaller. 

4. THE PROTOSOLAR AND THE PRESENT SOLAR VICINITY 

ABUNDANCES 

To compare GCE models of Y, C, and O with observations we need to use 
the best abundance determinations available for these elements. We consider 
that the two most accurate Galactic Y, C, and O determinations are the 
protosolar values (Asplund et al. 2009) and the M17 H n region values (Paper 
I)- 

4.1. Protosolar abundances 

In order to study the time-agreement with the protosolar values, we show 
in Figure 3 the predicted evolution of AY vs AO and AC vs AO from to 8.5 
Gyr (the Sun- formation time). The AY value for Figure 3 (and throughout 
this paper) is given by AY = Y — Y v , where Y v i s the primordial He abundance, 



and amounts to 0.2477 iPeimbert et al] (|2007l) . In this figure AO = O, and 



AO = O because at the time of the primordial nucleosynthesis O and C are 
not produced. For t = the models start at AY, AO, and AO equal to zero, 
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Fig. 3. AY vs AO (upper panel) and AC vs AO (lower panel) for the solar vicinity 
(r = 8 kpc). Evolution from t = (AY = AC = AO = 0) to 8.5 Gyr (Sun-formation 
time) predicted by models that as sume HWY (continuo us lines), and LWY (dashed 
lines) . : Protosolar values by lAsplund et al. (|2009h . Op en circles: Protosolar 
values from photospheric data by different author s compiled by Asplund et al. I |2009l . 
Table 4) and corrected for gravitational settling ( Asplund et al. 20091 ). 



and in this figure the evolution of the models stops at t = 8.5 Gyr, the time 
the Sun was formed. Also in this figure we present some of the most popular 
Y, C, and O protosolar abundances of the last 21 years. These photospheric 
solar abundances were obtained by Anders & Grevesse (1989), Grevesse & 
Noels (1993), Grevesse & Sauval (1998), Lodders (2003), Asplund, Grevesse 
& Sauval (2 005), and Lodd e rs, Pa lme & Gail (2009), abundances that were 
compiled by Asplund et al. ( 20091) in their Table 4. To obtain the protoso- 
lar abundances of the heavy elements, the photospheric abundances were in- 
creased by 0.04 dex, amount that takes into account the effect of gravitational 



settlin g. We used the protosolar X values shown in Table 4 by I Asplund et al 



(2009) to change the abundances by number to abundances by mass. 

From Figure 3 it can be noted that: a) the HWY model fits very well the 
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Y and O protosolar values, while the predicted C is about 1 a error higher 
than observed, b) the LWY model matches the Y and O protosolar values 
within 1 <t, but the predicted C is about 2 a lower than observed, c) the solar 
abundances predicted by our models are in much bett er agreement w i th the 



recent He, C, and O protosolar values determined by lAsplund et al 



(l2009h 



than with the previous ones compiled by them, and d) the He protosolar abun- 
dance determinations have remained almost constant over the years, while the 
C and O ones have decreased. 



4.2. Ay vs AO evolution compared with M17 and young B stars of the solar 

vicinity 

In Figure 4 we present the evolution of AY versus AO. The models 
presented in this figure for r equal to 6.75, 8, and 17 kpc are for the to 13 
Gyr range. The data should be compared with the end point of the evolution 
of the corresponding evolutionary track. In Figures 4, 5 and 7 the evolutionary 
tracks and the related observational data have the same color. 

Since there are no good abundance determinations for the outer parts of 
the Galactic disk, we test our models with data of irregular galaxies. There- 
fore, in Figure 4 we show the AY and AO values of the metal poor extragalac- 
tic H II regions determined bv lPeimbert et al.l (|2007l ) . Specifically, we compare 
our 17kpc track with NGC 2366, because it contains one of the brightest H n 
regions of the Galactic vicinity. The gaseous O abundances of the metal poor 
H II regions were increased by 0.10 dex to incl ude the fraction of O atom s 
embedded in dust grains inside the H n regions ( Peimbert fc Peimbert 2O10h . 



From that figure it can be noticed that the predicted chemical evolution of 
the Galaxy for large galactocentric radii (r > 17 kpc) at t = 13 Gyr behaves 
like irregular galaxies at present-time. 

In Fig ure 4 we also present th e abundances of B stars of the solar vicinity 
de rived by Przvbilla et al.l (120081 ) . We have adopted the Orion X value shown 



Carigi et al.l (|2006l ) for converting Xi/H (the abundance ratio by number 



of any i element) of B stars to Xi by mass. The AY abundance derived from 
B stars is in good agreement with the HWY and LWY models for r = 8 kpc. 
On the other hand the AO value is about 1.5 a and about 4.5 a smaller than 
th e values pred i cted b y the HWY and LWY models. The O/H values derived 
by Simon-Diaz ( 2010h for the B stars of t he Or ion star forming region are in 
good agreement with the iPrzvbilla et al.1 (|2008l ) determinations. 

We decided to include as observational constraints the O and He abun- 
dances of the H II region M17. This object has the best He abundance de- 
termination available because its degree of ionization is very high and the 
correction for neutral helium, that is always indirect, is the smallest for the 
well observed Galactic H II regions (see Paper I) . We show in Figure 4 two sets 
of AO values, one derived from O recombination lines, and the other derived 
from the O forbidden lines under the assumption that i 2 = 0.00. 

The point with the highest AY and AO values of the HWY model at 6.75 
kpc presented in Figure 4 (the point at t = 13 Gyr), should be compared 
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Fig. 4. Evolution of AY vs AO. The HWY model is shown in the upper panel and the 
LWY model in the lower panel. Chemical evolution tracks from (AY = AO = 0) 
to 13 Gyr (present time) for r = 6.75, 8, and 17 kpc (dotted-red, solid-green, and 
dashed-cyan lines, respectively). The three tracks in each panel partially overlap, 
from top to back they lie as follows: 17kpc, 8kpc, and 6.75kpc. The values for the 
H II regions and the B stars correspond to 13 Gyr. The M17 H II region at r = 
6.75 kpc is represented by a filled circle, for t 2 = 0.036 ( RL) and by an open ci rcle 



for t 2 = 0.00 (FL). The solar vicinity B stars value by IPrzvbilla et all l|2008l ) 



represented by a star. The extragal actic low metallicity H II regions are represented 
by filled triangles for t 2 0.00 (FL) jPeimbert et al.ll2007f ). A Y p = 0.2477 is adopted 
( Peimbert et al.ll20071 ). 



with the two M17 values. The AY model value is in good agreement with the 
two M17 values. The AO value for t 2 = 0.00 is more than 5a smaller than 
the model prediction, while the t 2 ^ 0.00 point is less than 2a higher than 
the AO model prediction. For the LWY model the predicted AY value is in 
good agreement with the two M17 values; but the AO value of the t 2 = 0.00 
point is 12cr smaller than the model prediction, while the AO value of the 
t 2 7^ 0.00 point is less than 2a smaller than the model prediction. Based on 
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this comparison we conclude that it is possible to get an excellent agreement 
with the t 2 ^ 0.00 AO value using a model with intermediate yields between 
the HWY and the LWY models, but that it is not possible to find a reasonable 
model to fit the AO value derived assuming i 2 = 0.00. 

Since the HWY and the LWY models fit the protosolar values and produce 
a reasonable fit to the observed M17 AO and AY values, we conclude that 
the protosolar values in the context of Galactic chemical evolution provide a 
strong consistency check to the O recombination abundances derived from H n 
regions and are in disagreement with the O abundances derived from forbidden 
lines under the assumption that t 2 = 0.00. Ad ditional support for this result 
is provided by ISimon-Diaz &: Stasinskal (]201lf ) who based on a study of 13 
B-type stars of the Orion star forming region OBI find that the stellar O/H 
abundances agree much better with the Orion H n region abundance derived 
from recombination lines than with the one derived from collisionally excited 
lines. 



4.3. AO vs AO evolution compared with H II regions, and young B, F, and 

G stars of the solar vicinity 

In what follows we will discuss the AO and AO observational values 
adopted from the literature and we will compare them with the HWY and 
LWY model predictions. In Figure 5 we show the 0-13Gyr evolution of the 
AO-AO relation for r = 7, 8, and 17 kpc. 

The O/H ratio of the LWY model for r = 17 kpc and t = 13 Gyr, cor- 
responds to the O/H value of NGC 2363, the best extragalactic metal-poor 
H II region for our purpose. This model is in good agreement with the C/O 
observed ratio. Similarly the O/H ratio of the HWY model for r = 17 kpc and 
t = 13 Gyr, corresponds to the O/H value of NGC 2363, this model is also in 
good agreement with the C/O observed ratio. The AO gas and dust values 
for NGC 2363 are those presented in Figure 4, while the AO value includes 
gas (jEsteban et al.l 120091) and dust compo nents. We have considered the Y 
and O/Z values from lPeimbert et al. ( 2007 ) to calculate X for NGC 2363 and 
to change the abundances by number to abundances by mass. 

For r ~ 8 kpc and t = 13 Gyr there are several reliable C and O abundance 
determinations of the solar vicinity. To compare our models with observations 
we chose: a) the average C/H and O/H values of NGC 3576 and Orion H II 
regions at r= 7.46 and 8.40 kpc and the correction for the fractions of C and 
O embedded in dust gra ins, b) the young F and G dwarf stars studied by 
Bensbv fc Feltzind J2006I) . and c) the average C and O values of B stars by 
Przvbilla et al.l |2008l ). 

In Figure 5 we compare our models with the average AO and AO values 
of NGC 3576 and Orion, and find a reasonable agreement. To convert the 
Xi /H values by numbe r to Xi by mass we adopted the X(Orion) value shown 
bv lCarigi etaH tod) . 



Bensbv fc FeltzinsJ (|2006l ) studied 51 F and G dwarf stars of the solar 



vicinity, 35 belonging to the thin disk and 16 to the thick disk. To compare 
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Fig. 5. AC vs AO diagram. The dotted red line represents the chemical evolution 
track from (AC = AO = 0) to 13 Gyr (present time) for r = 7 kpc, the other 
lines as in Figure 4. A description of the symbols follows. Filled pentagon: Average 
values of the M17 and M20 H II regions at r = 6.75 and 7.19 kpc, respectively. 
Empty sta r. Average values of youn g F-G dwarf stars of the solar vicinity, r ~ 8 
kpc, from bensbv fc Feltzin j (|2006h (see text). Filled square: Average values of 
NGC 3576 and Or ion H H regions at r = 7.46 and 8.40 kpc. Filled star: Average 
value of B stars bv lPrzvbilla et all (120081 ). Filled triangle: NGC 2363 H n region in 
a very metal-poor irregular galaxy (see text). 



the abundances predicted by our models with those of the youngest stars 
presumably recently formed, it is necessary to select the youngest subset of the 
thin disk stars. The metal richest stars are expected to be the youngest ones. 
We made three subsets of the metal richest stars of the thin disk containing 
4, 8 and 16 objects respectively, and obtained 12+log(0/H) average values 
of 8.87, 8.84, and 8.83 respectively, and 12+log(C/H) average values of 8.63, 
8.64, and 8.59 respectively. In Figure 5 we present the subset of 8 stars as 
repres entative of the pres ent day ISM, where we have adopted the A (Orion) by 
mass ( Carigi et al. 20061 ) for converting Xi/H to Ai by mass. The agreement 
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is very good. As we saw above the other two subsets of F and G stars produce 
similar C/H and O/H values also in good agreement with our models. 

To convert the C and O abun d ances by number to abundances by mass 
of the B stars by Przvbilla et al.l (|2008l) we have taken the Orion X value 



shown in ICarigi et al.l (|2006l ). From Figure 5 it can be seen that the AC 
value derived for the B stars is many sigma smaller than predicted by the 
HWY and the LWY models. Moreover it is also considerably smaller then 
the values derived from the F-G dwarf stars, and the average of NGC 3576 
and the Orion nebula. Maybe part of the reason for the smaller C abundance 
in the B stars can be du e to rotational mixing (e.g. Mevnet fc Maeder 2000t 



Fierro fc Georgievl 120081 ). The B stars AO value is about 1.5a and about 



4.5a smaller than the predicted HWY and LWY model values respectively. 
Furthermore the B stars AO value is in fair agreement with the F-G dwarf 
value and about 1.5er smaller than the the average value of NGC 3576 and 
Orion. These differences between the B stars and the other objects of the 
solar vicinity should be studied further. 

For r = 7 kpc and t = 13 Gyr we took the average of the abundances of 
the M17 and M20 H n regions as the o bservational constraint, the C/H and 
O/H gaseous values were obtained from lGarcia-Roias fe Estebanl (|2007j) and 
Paper I using recombination lines of C and O. We have adopted the M17 X 
value, obtained in Paper I, to convert the abundances by number of M17 and 
M20 to abundances by mass. Unfortunately there are no observational data 
representative of the chemical past at r = 7 kpc. 

The HWY model adjusts very well the current C and O values for r = 7 
kpc and marginally the present-day O for r = 8 kpc, but does not explain the 
current C value for r — 8 kpc. Alternatively the LWY model adjusts quite 
well both the C and O values for the solar vicinity, but predicts a higher O 
abundance for r = 7 kpc. 

In Figures 4 and 5 the evolutionary tracks partially overlap, for a given 
AO value the AY and AO values of the two different tracks are almost the 
same, but the time is very different. In Figure 6 we present the evolution of 
AO, AY, and AO as a function of time to appreciate the differences between 
the 7kpc and the 17kpc tracks. It is important to present the behavior of 
/i = Mg as /Mtot because it drives the O evolution. For example the 17kpc 
track reaches a AO value of about 2 x 10~ 3 at 13 Gyr, while the 7kpc track 
reaches this value at less than 1 Gyr, where the same AO value for both tracks 
occurs at a similar /i value. The small differences in AO, AY, and AO at 
a given /i value come from: a) the delay in the C and He enrichment of the 
ISM due to LIMS, b) the shape of the SFR, and c) the differences between 
the HWY and LWY. In Figure 6 we present the SFR behavior where it can 
be seen that for r = 7 kpc the SFR is not only different in shape to the one 
at r — 17 kpc but it is from one to two orders of magnitude higher. 

From the HWY model it has been found that about h alf of AY and AO 
have been formed by LIMS and half by MS, (see Figure 3 of Carigi et al.ll2005 , 
and Figure 6 of Paper I), while most of the AO has been formed by MS. LIMS 
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Fig. 6. Evolution of M gaa /Mtot = star formation rate, AO, AY, and AC for 
r = 7 kpc (dotted lines) and r = 17 kpc (dashed lines), assuming HWY (thick 
tracks) or LWY (thin tracks). Star formation history is in M Q pc~ 2 Gyr^ 1 units and 
the SFR for r = 17 kpc has been multiplied by 10. 



produce a v ery small amount of O acco rding; to yields by different authors (see 
Figure 7 of lKarakas fc Lattanzioll2007l ). in this paper we are using the yields 
by Marigo and collaborators. For the LWY model the fractions of AY" and 
AG due to MS are similar but smaller than for the HWY model. The similar 
fraction of C and He produced by MS and LIMS is also responsible for the 
behavior in Figures 4 and 5 where it can be seen that the AY/ AO and the 
AG/ AO relations are similar for the HWY model and for the LWY model for 
AO < 5 x 10" 3 . 



5. INTERMEDIATE WIND YIELDS 



From Figures 1 — 5 we conclude that the HWY and LWY models are in 
very good agreement with some of the data, but for other data the agreement 
is only fair. To improve the agreement with the observations we suggest the 
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Fig. 7. Chemical evolution model that assumes IWY=(HWY+LWY)/2.0. The left 
panels show the 0-13Gyr evolution of AY vs AO and of AC vs AO for r = 7 kpc 
(dotted red lines) and r — 8 kpc (thin solid green line), data as Figs. 4 and 5; the 
thick yellow lines show the evolution from to 8.5 Gyr for r = 8 kpc, data as in 
Fig. 3. The right panels show the present-day ISM abundance ratios as a function 
of galactocentric distance, data as in Fig. 1. 



use of intermediate wind yields (IWY) for the computation of the chemical 
evolution models. We define the IWY as the average of HWY and LWY. 

In Figure 7 we present the chemical abundance predictions derived from 
the IWY model for the most critical observations. In particular, it is remark- 
able the fit of the model with the observed values for: a) the present-day (t = 
13 Gyr) C/O gradient in the 6 < r(kpc)< 11 range with the gradient derived 
from Galactic H n regions, b) the t = 13 Gyr Ay, and AO values for r = 7 
kpc with those of M17, c) the t = 13 Gyr AC, and AO values for r = 7 kpc 
with the average of the M17 and M20 values, d) the t = 13 Gyr AO, and AO 
values for r = 8 kpc with the average values for NGC 3576, and Orion, and 
for young F and G stars of the solar vicinity, and e) the t — 8.5 Gyr AY, AO, 
and AO values for r = 8 kpc with the protosolar values. 
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TABLE 1 
EXTRAGALACTIC H II REGIONS 



Host Galaxy 


Type of Galaxy 


kj Dject 


v (kpc) 


1 0_l_l^rv/ C\ /W\ a 

lA+LOg(KJ/ tl) 


iz+log^/rij 




M 101 


Scdl 


H1013 


5.50 


8.85 ± 0.09 


8.67 ± 0.12 


-0.08± 0.15 






NGC 5461 


9.84 


8.61 ± 0.06 


8.30 ± 0.20 


-0.21± 0.22 






NGC 5447 


16.21 


8.64 ± 0.06 


8.20 ± 0.12 


-0.34± 0.14 






NGC 5471 


23.45 


8.35 ± 0.15 


7.79 ± 0.19 


-0.46± 0.24 


M 33 


ScdII-III 


NGC 595 


2.87 


8.81 ± 0.05 


8.63 ± 0.12 


-0.18± 0.13 






NGC 604 


4.11 


8.72 ± 0.03 


8.40 ± 0.11 


-0.22± 0.12 


M 31 


Sbl-II 


K932 


16.0 


8.74 ± 0.03 


8.49 ± 0.13 


-0.18± 0.14 


NGC 2403 


ScdIII 


VS44 


2.77 


8.73 ± 0.04 


8.32 ± 0.18 


-0.31± 0.19 



a Gaseous value plus dust correction. 



The fits presented in Figure 7 imply that the effects of migration have not 
been very important in the history of the chemical evolution of the Galaxy. A 
si milar conclusion on th e migration effects for the solar vicinity has been found 
bv lNavarro et al.l (|2010h based on the kinematic properties and the metallicity 
of thin disk stars. 



6. POSSIBLE IMPLICATIONS FOR OTHER SYSTEMS 

In is important to study how general is the chemical evolution model de- 
rived for the disk of the Galaxy and to explore its relationship with other 
spiral galaxies and the transition region between the disk and the bulge of 
the Galaxy. In what follows we present a preliminary discussion of these two 
topics. 



6.1. Other spiral galaxies 

Since the IWY model successfully reproduces the He, C, and O abundances 
in the Galactic disk and in the solar vicinity, we want to extend our study to 
other spiral galaxies with abundances determined by recombination lines in 
H II regions. 

In Table 1 we collect the gaseous abundances of 8 cxtragalactic H n re- 
gions that belong to four spiral galaxies: M31, M33, M101, and NGC 2403 
(jEsteban et al]|2002l [2009h . where we have included dust corrections identical 



to those considered for Galactic H n regions. Consequently, these data are 
consistent with the Galactic data used in our previous figures (e.g. Figure 7). 
In Table 1 , we show the name and type of the host galaxy, and the name and 
the galactocentric distance of each extragalactic H II region. 

In Figure 8 we show the present-day C/O and O/H values for r = 4,5, 15, 16 
kpc (from right to left) predicted by the IWY model of the Galactic disk and 
the extragalactic H n regions values shown in Table 1. 
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Fig. 8. C/O-O/H values at 13 Gyr predicted by IWY model (long-dashed lines). 
Small black circles: model predictions at different Galactocentric distances, from 4 
kpc (right) to 16 kpc (left). Observed abundances ratios of extragalactic H n regions 
in spiral galaxies: cyan squares: M101, big green circles: M33, blue triangle: M31, 
and magenta pentagon: NGC 2403. 



We notice that in the C/O-O/H relation presented in Figure 8: a) the 
spiral galaxy disks have a similar behavior to that of the Milky Way disk, b) 
all the extragalactic H n regions of the sample are well reproduced by IWY 
models, and c) we see the effect of the Z dependent yields for 12+log(0/H) 
> 8.4. More determinations of the O/H and C/H ratios in extragalactic H II 
regions with 12+log(0/H) values < 8.4 and > 8.8 are needed to test the yields 
further. 

Since the general C/O-O/H trend of the extragalactic H n regions can be 
addressed using the results of the GCE model, we can suggest that the four 
extragalactic spiral galaxies have a similar IMF, no selective outflows, and 
probably an inside-out formation scenario like our galaxy. 

It should be clear that Figure 8 does not correspond to models of the four 
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spiral galaxies, and that the r values correspond to our Galaxy and not to the 
other galaxies. 

In Table 2 we present some galactic information of M33, M101, and the 
Milky Way, the only galaxies with determined chemical gradients based on 
recombination lines. Specifically, we show the galactic photometric radius to 
25 mag per square second ( Rp) and the chemical gradients no rmalized to Rq 
corrected by d ust depletion festeban et al.lbocd l2005l 120091). The MW R a 
was taken from lde Vaucouleurs fc Pence! (|l978l ). 

In Figure 9 we show the O/H, C/O, and C/O values as a function of 
r normalized to Rq for the 8 extragalactic and the 8 Galactic H 11 regions 
presented in Figs. 8 and 7, respectively. Moreover we represent the chemical 
gradients shown in Table 2. 

From Figure 9, it can be noted that: i) the O/H slope for the three galaxies 
is the same within the errors, which suggests that for a given galaxy the 
O/H gradient normalized to Rq is a more meaningful representation than 
the standard gradient; ii) the O/H gradients extend from r/Ro ~ 0.2 to 
r/Ro ~ 1.0; iii) the O/H ratio at a given r/Ro differs by a constant among 
M33, M101 and the MW; and iv) for 12+log(0/H)~8.8 the slopes of the C/H 
and C/O gradients become steeper. 

The parallel behavior of the O /H gradients when they are plotted relative 
to r/Ro is remarkable, suggesting that the main mechanisms involved in the 
formation and evolution of spiral galaxies are similar. 

The slope of the O/H gradients is similar, but the O/H absolute value 
at a given r/Ro is different for each galaxy. In particular the MW shows 
a ratio 0.24 dex higher O/H than the average of M33 and M101 at a given 
t/Rq. Considering that the MW is an Sbc galaxy and M33 and M101 are Scd 
galaxies, this result probably implies that the earlier the type of a galaxy the 
higher the O/H value at a given r/Ro, this result might be due to the fact 
that the earlier the bulk of star formation the fainter the stellar luminosity 
at a given O/H value and consequently the higher the O/H value at a given 
i?o value. Considering that the MW is an Sbc galaxy and M33 and M101 are 
Scd galaxies, this result might be of a general nature and suggests that the 
earlier the type of a galaxy the higher the O/H value at a given r/Ro- More 
high accuracy determinations of O gradients in spiral galaxies are needed to 
test this result. 

From Figures 8 and 9 it follows that we need: i) to increase the r/Ro 
coverage, mainly in the shorter galactocentric distance, with abundance de- 
terminations of high accuracy, ii) to increase the sample by including spiral 
galaxies of different types, and iii) to make specific models for each galaxy. 
All in order to sort out possible differences in the galactic formation and evo- 
lution relative to that of the Galaxy an d to understan d the peculiar behavior 
of metal-rich stars and H 11 regions, see lCarigil (|2008T ). and to be able to test 
stellar yields for Z > Z & . 

The present paper only includes specific models for our galaxy, we present a 
preliminary discussion on the probable relevance of our results to the study of 
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TABLE 2 

ABUNDANCE GRADIENTS NORMALIZED TO R a a 



Galaxy 


i?o(k P c) b 


12+log(0/H) 


12+log(C/H) 


log(C/0) 


M101 


28.95 


-0.492 x(r/i? )+8.84 


-1.32x(r/i? )+8.99 


-0.65x(r/i? )+0.03 


M33 


6.83 


-0.492x(r/i? )+9.00 


-0.72x(r/i? )+8.93 


-0.22x(r/i? )-0.10 


MW 


11.25 


-0.495x(r/i? )+9.16 


-1.16x(r/i? )+9.50 


-0.65x(r/i? )+0.34 



1 -Ro = galactic photometric radius to 25 mag per square second. 
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Fig. 9. Chemical abundance ratios of H n regions normalized to galactic photometric 
radius (Ro) for each galaxy (see Tables 1 and 2). Milky Way: red asterisks (see Fig. 
7). Others spiral galaxies: symbols as Fig. 8. Lines represent the galactic disk 
gradients of M101, M33, and the Milky Way. 



other spiral galaxies. It is beyond the scope of this paper to produce models 
for other galaxies and therefore to present a detailed comp a rison with the 
results derived by other authors, e.g. IChiappini et al.l (|2003bl ); iMolla fc Diaz 
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Fig. 10. Evolution for r = 3 kpc predicted by the IWY model (long-dashed black 
lines). Star formation history normalized to 70 MqPC~ 2 Gyr -1 . C/O-O/H, O/Fe- 
Fe/H, and time-Fe/H relations from to 7 Gyr. Area enclosed by dot- short- dashed- 
green lines: Galactic bulge red giant stars presented by Cescutti et al. (2009). 
Filled symbol s : micr olensed dwarf and subgiant stars in the Galactic bulge from 
Bensbv et al.l (|2010al ). red squares: old stars with mean age of 9.8 ± 4.0 Gyr, blue 
triangles: young stars with mean age of 4.5 ± 1.6 Gyr. Ages were normalized to 



the ag e of the models (13 Gyr.) Dotted lines: protosolar values from lAsplund et al 
1 20091) . 



( 20051); iRenda et al 



|2005l) ; lYin et al. I (|2009h : iMarcon-Uchida etTaO |2010h : 
Magrini et al.l ()2010f ) . using different sets of data. In future pape rs we plan 
to present detailed models for M31 (jMeneses-Govtia et al.1 120111 ) and M33 
(Robles-Valdez, Carigi, & Bruzual 2011 in preparation). 



6.2. The transition region between the inner disk and the Galactic bulge 



A fraction of the stars observed in the direction of the Galactic center 
might have been formed in the true bulge and another fraction in the inner 
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Galactic disk. Therefore we considered interesting to compare our IWY model 
for the innermost part of the disc with the abundances of these stars. 

Since the inside-out model adopted as the Galactic formation scenario in 
this paper diverges for r = 2 kpc (see Section 2), we obtained the chemical 
evolution of the innermost ring of the Galactic disk that our GCE model is 
able to consider, r = 3 ± 0.5 kpc (with a 1 kpc width). The Galaxy at this 
radius formed efficiently: during the first Gyr the halo formed with a collapse 
time-scale of 0.5 Gyr and then the disk formed with a formation time scale of 
1.0 Gyr. At the end of the evolution (t = 13 Gyr) the gaseous mass amounts 
to 7% of the baryonic mass inside the r = 3 ± 0.5 kpc ring. 

In Figure 10 we show the C/O-O/H, O/Fe-Fe/H, and Fe/H evolution ob- 
tained with the IWY model between and 7 Gyr. We include the evolution 
from to 7 Gyr and exclude the 7-13 Gyr range because at 7 Gyr the model 
reaches the highest Fe/H value of the stellar data. Moreover, in the excluded 
range the SFR is very low and the probability to form stars younger than 6 
Gyr (in the 7 to 13 Gyr range) is negligible. 

In Figure 10 (C/O-O /H panel) we show also the area covered by the data 



of ICescutti et al.l (|2009l ). These authors collect ed [C/O] and fO/Hl r atios 



for the Galactic bulg e red giants determined by iFulbright et al.l (|2007h and 



Melendez et al.l (|2008l ). and used the observed data to obtain t he initial stellar 



C/H a nd O/H values relative to the solar abundances due to lAsplund et al 



(|2005l ). In order to convert the Cescutti et al. data to abundances by numbe r 
we used, only in this figure, the solar abundances bv lAsplund et al. (2005 s ) 



For th e other figures, we assumed the protosolar abundances bv lAsplund et al 
|2009h . 

A description of the model in the C/O-O/H panel of Figure 10 follows. 
During the halo formation, the IWY model predict an increase of 12+log(0/H) 
to 8.6 dex and an increase of 12+log(C/0) to -0.7 dex . After the halo for- 
mation finishes, at t = 1 Gyr, the O/H ratio decreases due to the dilution 
produced by the enormous amount of primordial gas that is accreted to form 
the disk. Later on this accretion causes a rapid SFR increase, producing an 
O/H increase. The C/O rise from -0.7 dex to -0.5 dex is due mainly to LIMS, 
while the C/O rise from -0.5 dex to 0.0 dex is due to both, massive stars and 
LIMS. 

In Figure 10 (O/Fe-Fe/ H and time-Fe/H panel s), we included the chemical 
abundances determined by iBensbv et al.1 (|2010al ) in microlcnsed dwarf and 
subgiant stars of the Galactic bulge. In order to convert the data by Bensby 
et al. to abundances by number, we used their solar abundances. The stellar 
ages were normalized to the age of the model (t = 13.0 Gyr). 

Our mode l at r = 3 kpc successfully reproduces the O/Fe-Fe/H obtained by 



Bensbv et alJ (|2010al L but the timc-Fc/H relation is only partially reproduced. 



Our model explains all old-metal-poor stars and two young stars, but not 
the oldest metal-rich star in their sample (see the square at 12+log(F e/H)^ 
7.9) which behaves like a common bulge star, see iBallero et al.l (|2007l ). The 
youngest and metal-rich stars would be disk stars that formed at r < 3 kpc. 
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Since Fe is produced by massive stars and binary systems of LIMS, to a first 
approximat ion Fe behaves like C, t hat is produced by MS and single LIMS, see 
for example I Akerman et al.l (|2004l ). the C/O-O/H discussion presented above 
can be used to explain the O/Fe-Fe/H and time-Fe/H panels of Figure 10. 



Bensbv et al.l (|2010bl ) have studied red giant stars in the inner Galactic 



disc and find that the abundance trends of the inner disc agree very well with 
those of the nearby thick disc, and also with those of the Galactic bulge. 
Based on the stellar results of the thick and thin dis ks of the solar vicinity, of 
the inner Galactic d isc, and of the Galactic bulge ( Bensbv fe Feltzing 20061 : 
Bensbv et al. 2010al fb). they sugge st that the bulge and t he disk have had 



similar chemical histories. Moreover Alves-Brito et al. I (l2010h also suggest that 
the bulge and local thick disk stars experienced similar formation timescales, 
star formation rates and initial mass functions. 
_ Cescutti et al. ( 20091 ) present a Galactic bulge model, computed by Ballero et al 



(|2007l ). that, in order to reproduce the total stellar mass, the iron distribution 
function, and the a/Fe-Fe/H relations (all constraints obtained from giant 
stars) they assume: a formation time scale of ~ 0.1 Gyr, as well as a star 
formation efficiency one order of magnitude higher and an IMF flatter for 
77i > 1 M Q than those considered in our models. This model, with a SFR as 
lo ng as 0.5 Gyr, cannot explain the ages of the youngest bulge dwarfs found 
bv lBensbv et all |2010al) . 

Note that, if we focus only on C/O-O/H and O/Fe-Fe/H relations deter- 
mined in dwarf stars of the Galactic bulge we cannot distinguish between a 
scenario with a very rapid infall, effi cient star form a tion, and a high relative 
formation of massive stars, that by ICescutti et al.l (|2009l ). and another sce- 
nario with an order of magnitude less rapid infall, an order of magnitude less 
efficient star formation, and a lower relative formation of massive stars (the 
one presented in this paper). 

Nevertheless, if we focus on the age-Fe/H relation shown by bulge dwarf 
stars, our model, with an extended bulge formation time and therefore an 
ex tended star formation history, produces a better fit to the data than that 
bv lCescutti et al.l (|2009h . 

In Figure 11 we show the iron distribution function predicted by the IWY 
model for r = 3 kpc in the Galactic disk from to 7 Gyr. Since there are no 
good abundance determinations for the inner parts of the Galactic disk, we 
compare our model with the 104 stars that belong to the outer most field along 
the bu lge minor axis at b=-12, see lower panel of Figure 14 bv lZoccali et al. 
|2008h . We have used 12+log(Fe/H) = 7.50 to convert the [Fe/H] values of 
the Zoccali et al. sample to Fe/H abundances by number. From that figure, 
we can say that a non-negligible fraction of the stars in the direction of the 
bulge might belong to the inner disk. Moreover, we could comment on that 
bulge sample: i) for 12+log(Fe/H)< 7 the contamination of the thick disk 
is high, ii) for 7 < 12+log(Fe/H)< 8 the contamination of the thin disk is 
lower and the true bulge stars dominate, and hi) for 12+log(Fc/H)> 8 the 
contribution of the metal- rich stars of the thin disk is important. 
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Fig. 11. Iron distribution function. Predictions for r = 3 kpc model that assume 
IWY (long-dashed black line) until 7 Gyr. Observed distribution for the outermost 
fields (b=-12) along the bulge minor axis bv lZoccali et al. (2008) (continuous green 
line). Dotted line: protosolar Fe/H value from Asplund et al. ( 20091 ) . 



Note that in Figure 11 we are using only the outermost 104 stars of the 
total 800 stars in the dir ection of the bulge, and that our model is for r = 3 



kpc. On the other hand Cescutti fe Matteuccil (j 2 1 lh show that their bulge 



model (the same as Cescutti et al. ( 20091 )) follows the constraint provided by 



the 800 stars very well using a Salpeter IMF. The model by Cescutti and 
Matteucci predicts a smaller number of stars for [Fe/H] values higher than 
solar, probably indicating that the bulge sample has been contaminated by 
innermost disk stars. 

From the previous discussion it follows tha t the Galactic bulge is a complex 
structure that should be studied further (e.g. Zoccali 20lot) . 



7. CONCLUSIONS 



We have made models with three different sets of yields that differ only 
on their Z dependence at solar metallicities for massive stars: HWY, IWY, 
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and LWY. The HWY and the LWY have been used before by us, while in 
this paper we introduce the IWY, that are given by (HWY + LWY)/2. We 
find that the IWY galactic chemical evolution models produce better fits to 
the observational data than either the HWY or the LWY galactic chemical 
evolution models. 

We present a Galactic chemical evolution model based on the IWY for the 
disk of the Galaxy that is able to fit: a) the C/O vs O/H, C/Fe vs Fe/H, O/Fc 
vs Fe/H, and Fe/H vs t relations derived from halo and disk stars of different 
ages in the solar vicinity, b) the O/H, C/H, and C/O abundance gradients 
(slopes and absolute values) derived from Galactic H II regions, c) the He/H, 
C/H, O/H, Fe/H protosolar abundances, and d) the He/H and O/H values of 
the galactic H n region M17. 

We find that in general about half of the freshly made helium is produced 
by massive stars and half by LIMS, and that a similar situation prevails for 
carbon, while most of the oxygen is produced by massive stars. The agreement 
of the He/O and C/O ratios between the model and the protosolar abundances 
implies that the Sun formed from a well mixed ISM. 

We note that the agreement of our model with the protosolar abundances 
and the Sun-formation time supports the idea that the Sun originated at a 
galactocentric distance similar to that of the solar vicinity. 

We show that chemical evolution models for the Galactic disk are able to 
reproduce the observed AY and AO protosolar values and the AY and AO 
values derived for Ml 7 based on H, He and O recombination lines, but not 
the M17 AY and AO values derived from T(4363/5007) and O collisionally 
excited lines under the assumption of t 2 = 0.00. This result provides a con- 
sistency check in favor of the presence of large temperature variations in H II 
regions and on the method based on the H, He, C and O recombination lines 
to derive abundances in H n regions. 

We obtain that the IWY chemical evolution model of the Galactic disk 
for the present time, in the galactocentric range 6 < r(kpc) < 11, produces 
a reasonable fit to the O/H vs C/O relationship derived from H II regions 
of nearby spiral galaxies. The yields predict an increase of the C/O ratio 
with O/H starting from 12+log(0/H)^8.4 that is observed in our Galaxy 
and in nearby galaxies. The O/H vs C/O relationship might imply that spiral 
galaxies have a similar IMF, no selective outflows, and probably a formation 
scenario similar to that of our galaxy. 

We find a remarkable parallelism of the O/H gradients for M33, M101, and 
the Galaxy when they are plotted with respect to r/Ro (Ro is the galactic 
photometric radius to 25 mag per square second), suggesting some common 
mechanisms in the formation and evolution of spiral galaxies. The O/H ratio 
at a given r/Ro differs by a constant among M33, M101 and the MW. The 
MW shows a 0.24 dex higher O/H ratio than the average of M33 and M101 
at a given r/Rg. 

We also find that the results for our model at r — 3 kpc can explain: 
a) the C/O-O/H and O/Fe-Fe/H, relations, and b) partially the Fe/H-time 
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relation and the Fe distribution function derived from stellar observations in 
the direction of the Galactic bulge. We find that stars belonging to the thin 
and thick discs make a significant contribution to these relations. 

Future work to advance in this subject requires: a) to advance in the 
study of the galactic bulge to be able to quantify the stellar contributions 
due to the inner disk and the true bulge, b) to increase the H n regions r/Ro 
coverage, mainly in the shorter galactocentric distance, with high accuracy 
C/H and O/H abundance determinations; c) to increase the sample of spiral 
galaxies of different types with O /H gradients of high accuracy, and d) to make 
specific models for each galaxy. All in order to sort out possible differences in 
the galactic formation and evolution of other galaxies relative to that of the 
Galaxy. 

We thank Jorge Garcia- Roj as and Cesar Esteban for useful discussions. 
We are also grateful to the anonymous referee for a careful reading of the 
manuscript and many excellent suggestions. This work was partly supported 
by the CONACyT grants 60354 and 129753. 
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